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Abstract. The expression for the density operator of the damped harmonic oscillator is derived from 
the master equation in the framework of the Lindblad theory for open quantum systems. Then the 
von Neumann entropy and effective temperature of the system are obtained. The entropy for a state 
characterized by a Wigner distribution function which is Gaussian in form is found to depend only on 
the variance of the distribution function. 



1 Introduction 



In the last two decades, the problem of dissipation in quantum mechanics, i.e. the consistent 
description of open quantum systems, was investigated by various authors [1-5] (for a recent 
review see ref. [6]). It is commonly understood [4, 7] that dissipation in an open system results 
from microscopic reversible interactions between the observable system and the environment. 
Because dissipative processes imply irreversibility and, therefore, a preferred direction in time, 
it is generally thought that quantum dynamical semigroups are the basic tools to introduce 
dissipation in quantum mechanics. In the Markov approximation the most general form of the 
generators of such semigroups was given by Lindblad [8]. This formalism has been studied for 
the case of damped harmonic oscillators [9, 10] and applied to various physical phenomena, for 
instance, the damping of collective modes in deep inelastic collisions in nuclear physics [11]. 
In [12] the Lindblad master equation for the harmonic oscillator was transformed into Fokker- 
Planck equations for quasiprobability distributions and a comparative study was made for the 
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Glauber P, antinormal ordering Q and Wigner W representations. In [13] the density matrix for 
the coherent state representation and the Wigner distribution function subject to different types 
of initial conditions were obtained for the damped harmonic oscillator. From a preceding analysis 
[14] resulted an analogy in form of the Lindblad master and Fokker-Planck equations with the 
corresponding equations from quantum optics, which is based on the quantum mechanics of 
the damped harmonic oscillator [4] and the corresponding Brownian motion. The equation of 
motion generally used in the theory of Brownian motion is the master equation for the density 
operator or the Fokker-Planck equation satisfied by the distribution function [4, 15, 16]. 

In the present study we are also concerned with the observable system of a harmonic oscillator 
which interacts with an environment. The aim of this work is to explore further the physical 
aspects of the Fokker-Planck equation which is the c-number equivalent equation to the master 
equation. Generally the master equation for the density operator gains considerably in clarity if 
it is represented in terms of the Wigner distribution function which satisfies the Fokker-Planck 
equation. Within the Lindblad theory for open quantum systems, we will describe the evolution 
of the considered system towards a final equilibrium state. First we derive a closed form of the 
density operator satisfying the Lindblad master equation and then calculate the von Neumann 
entropy and the effective temperature of the quantum-mechanical system in a state characterized 
by a Wigner distribution function which is Gaussian in form. The two quantities are shown to 
evolve to their equilibrium values. 

Physically, entropy can be interpreted as a measure of the lack of knowledge (disorder) of the 
system. The effective temperature associated with the Bose distribution has been introduced [17] 
in connection with the entropy obtained in the quantum theory of relaxation for the harmonic 
oscillator. The entropy for an infinite coupled harmonic oscillator chain has also been calculated 
for classical [18] and quantum mechanical systems [19] represented by a phase space distribution 
function. In the present work we extract the density operator of the damped harmonic oscillator 
in the Lindblad theory by using a technique analogous to those applied in the description of 
quantum relaxation [17, 19, 20, 21]. Denoting by pit) the density operator of the damped 
harmonic oscillator in the Schrodinger picture, the von Neumann entropy S(t) is given by the 
expectation value of the logarithmic operator In p : 

S(t) = —k < In p >= — /cTr(plnp), (1.1) 

where k is Boltzmann's constant. Accordingly, the calculation of the entropy reduces to the 
problem of finding the explicit form of the density operator. 

The content of this paper is arranged as follows. In Sec. 2 we write the master equation 
for the density operator of the harmonic oscillator. In Sec. 3 we derive a closed form of the 
density operator satisfying the master equation based on the Lindblad dynamics. By using the 
explicit form of the density operator, in Sec. 4 we calculate the von Neumann entropy and 
time-dependent temperature and analyze their temporal behaviour. Finally, concluding remarks 
are given in Sec. 5. 
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2 Master equation for the damped quantum harmonic 
oscillator 



L(p(0), (2.2) 



The rigorous description of the dissipation in a quantum mechanical system is based on the 
quantum dynamical semigroups [2, 3, 8]. According to the axiomatic theory of Lindblad [8], the 
usual von Neumann-Liouville equation ruling the time evolution of closed quantum systems is 
replaced in the case of open systems by the following equation for the density operator p: 

Here, $ t denotes the dynamical semigroup describing the irreversible time evolution of the open 
system in the Schrodinger representation and L the infinitesimal generator of the dynamical 
semigroup $ t . Using the structural theorem of Lindblad [8], which gives the most general form 
of the bounded, completely dissipative Liouville operator L, we obtain the explicit form of the 
most general time-homogeneous quantum mechanical Markovian master equation: 

dp(t) 
dt 
where 

L( P (t)) = ~[H, P (t)} + ^£(Kp(0,^] + Wi,Pit)V}]). (2.3) 

3 

Here H is the Hamiltonian of the system. The operators Vj and Vj are bounded operators on 
the Hilbert space H of the Hamiltonian. 

We should like to mention that the Markovian master equations found in the literature are 
of this form after some rearrangement of terms, even for unbounded Liouville operators. In this 
connection we assume that the general form of the master equation given by (2.2), (2.3) is also 
valid for unbounded Liouville operators. 

In this paper we impose a simple condition to the operators H, Vj, Vj that they are functions 
of the basic observables q and p of the one-dimensional quantum mechanical system (with [q, p] = 
ihl, where I is the identity operator on 7i) of such kind that the obtained model is exactly 
solvable. A precise version for this last condition is that linear spaces spanned by first degree 
(respectively second degree) noncommutative polynomials in q and p are invariant to the action 
of the completely dissipative mapping L. This condition implies [9] that Vj are at most first 
degree polynomials in q and p and H is at most a second degree polynomial in q and p. Then 
the harmonic oscillator Hamiltonian H is chosen of the general form 

H = H + ^(qp+pq), H = ^ + ^fq\ (2.4) 

With these choices the Markovian master equation can be written [6, 10]: 

dp % % % 

jj; = -ft[ H o,p\ ~ ^"(A + /i)[g,pp + Pp] + ^(A - p)[p,pq + qp] 

DpP -[l %P\\ - [P,P]] + ^([q, [P,P]} + [P, [?>]]), (2.5) 
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where D pp , D qq and D pq are the diffusion coefficients and A the friction constant. They satisfy 
the following fundamental constraints [6, 10]: 

i) D pp > 0, ii) D qq > 0, iii) D pp D qq - D 2 pq > \ 2 h 2 /A. (2.6) 

In the particular case when the asymptotic state is a Gibbs state 

p G (oo) = e-^/Tre-^, (2.7) 

these coefficients reduce to 

A + /i Tiuj A — /i h hcu 

D V v = fimuj coth — — , D aa = coth — — , D vo = 0, (2.8) 

pp 2 2kT' qq 2 mw 2kT J y J 

where T is the temperature of the thermal bath and the fundamental constraints are satisfied 
only if A > 

3 Density operator of the damped harmonic oscillator 

By introducing the real variables x±,x 2 corresponding to the operators q,p: 

Imuj 1 

I1 = W 12 = 73^' (3 ' 1) 

in [12, 13] we have transformed the master equation for the density operator into the following 
Fokker-Planck equation satisfied by the Wigner distribution function W(xi, x 2 , t) : 



9t i,,ti,2 sx t 2.^ * j ai ( ai 4 



where 



A-/i -w \ n w = \( "^A^ D pq \ 
uj X + I2j> W h\ D pq D pp /muj)- ^ 

Since the drift coefficients are linear in the variables X\ and x 2 and the diffusion coefficients 
are constant with respect to X\ and x 2 , Eq. (3.2) describes an Ornstein-Uhlenbeck process 
[22]. Following the method developed by Wang and Uhlenbeck [22], we solved in ref. [13] this 
Fokker-Planck equation, subject to either the wave-packet type or the 5-function type of initial 
conditions. 

In the present paper we consider the underdamped case < uj) of the harmonic oscillator 
[6, 10]. If the initial condition for the Fokker-Planck equation is of a Gaussian (wave-packet) 
type (xio and x 2 $ are the initial values of x\ and x 2 at t = 0, respectively) 

W w { Xl , x 2 , 0) = exp{-2[(x 1 - x 10 ) 2 + (x 2 - x 20 ) 2 ]}, (3.4) 
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the solution of Eq. (3.2) is given by [13]: 

W w (xi, x 2 ) = ^=== exp{--J- y> w (xi - xif + tp w {x 2 - x 2 f + Xw{xi - x 1 )(x 2 - x 2 )]}(3.5) 

7ThUJy/—B w B w 



where 



B w = 9l g 2 - -gl 9l = g* 2 = ^ + d j-(e^ - 1), g 3 = 2[e~™ + f (1 - e" 2 *)], (3.6) 

K = gia* 2 + 92a 2 ~ 9s, i>w = 9i + 92 - 9s, Xw = 2(#ia* + g 2 a) - g 3 (a + a*). (3.7) 

We have put a — (/i — A = —A — id, d\ = (a 2 mujD qq + 2aD pq + D pp /mw)/h, d 2 = 

(mwDgq + 2/iDpq/uj + D pp /muj) /h and fi 2 = oj 2 — p 2 . The functions X\ and x 2 , which are also 
oscillating functions, are found to be [13]: 

x 1 = e~ xt [(x 10 (cos fit + ^ sin fit) + ^20^ sin fit], (3.8) 

x 2 = e~ A *[(x 2 o(cosfit — ^ sin fit) — ^10^ sin fit]. (3.9) 

The solution (3.5) of the Fokker-Planck equation (3.2), subject to the wave-packet type of initial 
condition (3.4) can be written in terms of the coordinate and momentum (< A >= Tr(pA) 
denotes the expectation value of an operator A) as [6, 10]: 

W ^ q ' P) = 2^7! eXP ^~27^^ _ < <? >) 2 + ^(P- < P >) 2 - 2 x{q~ < q >)ip~ < V >)]},(3.10) 



where 



< q >= e ~ A '[(cosfit + ^ sin fit) < q(0) > sin fit < p(0) >], (3.11) 

fi mil 

2 

<p>= e - xt [ — sin fit < q(0) > +(cosfit - g sin fit) < p(0) >], (3.12) 

, „ 2 „ 2 hmu 3 . . 

<f) = a pp =<p >-<p>= — ^"Ptc, (3- 13 ) 

^ = a qq =<q 2 >-<q> 2 =-^^^ w , (3.14) 

X = o-pg = ^ < + > - < q ><P >= ^XtD, <^ = 0^-X 2 - (3-15) 

To obtain the explicit form of the density operator, we apply, like in [19, 21], the relation 
p = 2ithN {W s (q,p)}, where W s is the Wigner distribution function in the form of standard 
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rule of association and N is the normal ordering operator [23, 24] which acting on the function 
W s (q, p) moves p to the right of q. By the standard rule of association we mean the correspondence 
puiqn _^ ynpm k e t weeri functions of two classical variables (q,p) and functions of two quantum 
canonical operators (q,p). The Wigner distribution function (3.10), which is in the form of the 
Weyl rule of association [25] , can be transformed into the form of standard rule of association 
[26] by performing the operation: 

W s (q,p) = exp(^ih^)W(q,p). (3.16) 

We get the following Wigner distribution function: 

W s (q,p) = ^=exp{-^(q- < q >) 2 + ij(p- < p >) 2 (3.17) 

-2(X - i\){q~ <q>)(p-<p >)]}, (3.18) 

where 

^ = H-(x-^\) 2 ■ (3.19) 

The normal ordering operator N can be applied upon the Wigner function W s in Gaussian form 
by using McCoy's theorem [23, 24]: 

N[exp(Aq 2 + Bp 2 + Gqp)} = [Jexp(-^ 7 )] 1/2 exp(ag 2 + (3p 2 + 7 gp), (3.20) 

where a = A/C, (3 = B/C, C = sinhT/rj, T = -^( 7 2 - Aaf3) 1 ' 2 , with J = coshT + 
ihj sinhT /r = 1/(1 — ihG). After a straightforward, but lengthy calculation, we obtain the 
following expression of the density operator: 

p = - 7 =exp{-ln^- F -^- 7 =cosh (1 + ^) 

x[<f>(q- < q >f + ^{P- <P >) 2 -x[2(g- < q >)(p- <p>)-ih]]}. (3.21) 

This form is analogous to those obtained in [19, 21]. In particular, if D qq = D pq = and /i = A, 
we obtain the Jang's model [21, 27] on nuclear dynamics based on the second RPA at finite 
temperature. The density operator (3.20) has a Gaussian form, as expected from the initial 
form of the Wigner distribution function. While the Wigner distribution is expressed in terms 
of real variables q and p, the density operator is a function of operators q and p. When time 
t — > oo, the density operator tends to 

p(oo) = exp{-— \n 2 ^ + h {a pp (oo)q 2 + o qq (oo)p 2 - a pq (oo)(qp + pq)}}, (3.22) 



where a = a pp (oo)a qq (oo) — a pq (oo) and [6, 10]: 

(T " (o ° ) = 2mVA(W^) [mV(2A(A + /i) + Uj2)Dqq + Uj2Dpp + 2m ^ 2(A + (3 ' 23) 
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arm OO = 



' [m 2 uj i D qq +(2\(\- p)+oo 2 )D pp -2muo 2 (\- p)D pq l (3.24) 



2A(A- • <>-) 



<W°°) = 2mA(A 2 + ^ 2 ) [ ~ mV(A + ^ w + (A ~ + 2m(A ' ~ ^ 2)jDm] - (3 ' 25) 

In the particular case (2.8), 

°"oo(oo) = coth — — , am, (oo) = coth — — , <r BO (oo) =0 (3.26) 

and the asymptotic state is a Gibbs state (2.7): 

p G (oo) = 2smh — exp{-^(-p 2 + — , 2 )}. (3.27) 



4 Von Neumann entropy and effective temperature 

By using the relations (3.13-15), we get the expectation value of hip: 

< hip >= Ink - -\n(S - — ) - — In — = . (4.1) 

Denoting Uv = \/5 — h/2, we finally obtain the following expression of the von Neumann entropy: 

S(t) = k[(v+ l)hM> + 1) - v\nv\. (4.2) 

h 2 u 2 

Since 5 = — tpw^w, the function v becomes 



V = ^- B '-2- (4 - 3) 



where 



^ — , , . s,_w_ diet* f2 Mil c?9 „doi 

= exp(-4A t )(25Re^ r - _ + ^ - ^ + 2^) 

-2ex P (-2A i )[Re ((^ + ^)exp2« f ) - § + *] + M - fS. (4.4) 

It is worth noting that the entropy depends only upon the variance of the Wigner distribution. 
When time t -> oo, v tends to s = uj(d 2 2 /\ 2 - K| 2 /(A 2 + Vt 2 )) 1 / 2 /2Vt - 1/2 and the entropy 
relaxes to its equilibrium value S'(oo) = k[(s + 1) ln(s + 1) — s In s]. 

The expression (4.2) is analogous to those previously obtained [17, 18] in the theory of 
quantum oscillator relaxation and for the description of a system of collective RPA phonons 
[21]. It should also be noted that the expression (4.2) has the same form as the entropy of a 
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system of harmonic oscillators in thermal equilibrium. In the later case v represents, of course, 
the average of the number operator [19]. Eq. (4.2) together with the function v defined by (4.3), 
(4.4) is the desired entropy for the system. Although the expression (4.2) for the entropy has 
a well-known form, the function v induces a specific behaviour of the entropy. It is clear that 
the time dependence of the entropy is given by the damping factors exp(— 4At), exp(— 2Xt) and 
the oscillating function exp 2iQt. The complex oscillating factor exp 2iQt reduces to a function 
of the frequency cu, namely exp2iu;t for \i — > or if fi/Q <C 1 (i.e. the frequency uj is very large 
as compared to //). 

In the case of a thermal bath (2.7), (2.8), a time-dependent effective temperature T e can be 
defined [17, 21], by noticing that when t — > oo, v tends, according to (2.24) in ref. [13], to the 
average thermal phonon number < n >= (exp(hu/kT) — l) -1 . Thus v can be considered as 
giving the time evolution of the thermal phonon number, so that we can put in this case 

(exp^-1)- 1 ^. (4.5) 

The function v vanishes at t — 0. From (4.5) the effective temperature T e can be expressed as 

Accordingly, we can say that at time t the system is in thermal equilibrium at temperature T e . 
In terms of the effective temperature, the von Neumann entropy takes the form 

S(t) = ^ Hn[l-exp(--^)]. (4.7) 

T e(exp — -1) 

As t increases, the effective temperature approaches thermal equilibrium with the bath, T e — > T. 



5 Concluding remarks 



Recently we assist to a revival of interest in quantum brownian motion as a paradigm of quantum 
open systems. There are many motivations. The possibility of preparing systems in macroscopic 
quantum states led to the problems of dissipation in tunneling and of loss of quantum coher- 
ence (decoherence). These problems are intimately related to the issue of quantum-to-classical 
transition. All of them point the necessity of a better understanding of open quantum systems 
and require the extension of the model of quantum brownian motion. Our results allow such 
extensions. The Lindblad theory provides a selfconsistent treatment of damping as a possible 
extension of quantum mechanics to open systems. In the present paper we have studied the 
one-dimensional harmonic oscillator with dissipation within the framework of this theory. We 
have first obtained the explicit form of the density operator from the master and Fokker-Planck 
equations. The density operator in a Gaussian form is a function of the position and momentum 
operators in addition to several time dependent factors. Then the density operator has been 
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used to calculate the von Neumann entropy and the effective temperature. The temporal be- 
haviour of these quantities shows how they approach their equilibrium values. In a future work 
we plan to discuss the von Neumann entropy in association with uncertainty, decoherence and 
correlations of the system with its environment [28, 29]. 
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